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Abstract 



The modern status of basic low energy QCD parameters is reviewed. It is demonstrated, that 
the recent data allows one to determine the light quark mass ratios with an accuracy 10-15%. 
The general analysis of vacuum condensates in QCD is presented, including those induced by 
external fields. The QCD coupling constant a s {rn^) is found from the r-lepton hadronic decay 
rate. The contour improved perturbation theory includes the terms up to aj.. The influence of 
instantons on a s (m 2 ) determination is estimated. V-A spectral functions of r-decay are used for 
construction of the V-A polarization operator Hv-a(s) in the complex s-plane. The operator 
product expansion (OPE) is used up to dimension D=10 and the sum rules along the rays in the 
complex s-plane are constructed. This makes it possible to separate the contributions of operators 
of different dimensions. The best values of quark condensate and a s (0\qq\0} 2 are found. The 
value of quark condensate is confirmed by considering the sum rules for baryon masses. Gluon 
condensate is found in four ways: by considering of V+A polarization operator based on the r- 
decay data, by studying the sum rules for polarization operators momenta in charmonia in the 
vector, pseudoscalar and axial channels. All of these determinations are in agreement and result 
in {(a> s /ir)G 2 ) = 0.005 ±0.004 GeV A . Valence quark distributions in proton are calculated in QCD 
using the OPE in proton current virtuality. The quark distributions agree with those found from 
the deep inelastic scattering data. The same value of gluon condensate is favoured. 
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1 Introduction 



Nowadays, it is reliably established that the true (microscopic) theory of strong interaction is quantum 
chromodynamics (QCD), the nonabelian gauge theory of interacting quarks and gluons. The main 
confirmation of QCD comes from considering the processes at high energies and high momentum trans- 
fers, where, because of asymptotic freedom, the high precision of theoretical calculation is achieved 
and comparison with experiment confirms QCD with a very good accuracy. In the domain of low 
energies and momentum transfers (by such a domain in this paper I mean the domain of momentum 
transfers Q 2 ~ 1 — 5 GeV 2 ) the situation is more complicated: the QCD coupling constant a s is large, 
a s ~ 0.5 — 0.3 and many loops perturbative calculations are necessary. Unlike quantum electrodynamics 
(QED), the vacuum in QCD has a nontrivial structure: due to nonperturbative effects, non-zero fluc- 
tuations of gluonic and quark fields persist in QCD vacuum. The nontrivial vacuum structure of QCD 
manifests itself in the presence of vacuum condensates, analogous to those in condense matter physics 
(for instance, spontaneous magnetization). Therefore, a s corrections and nonperturbative effects must 
be correctly accounted in QCD calculation in this domain. 

At lower energies and Q 2 < 1 GeV 2 analytical QCD calculations are not reliable. The useful methods 
are: the chiral effective theory, lattice calculations and various model approaches. It is, however, very 
desirable to have a matching of all these approaches with QCD calculations at Q 2 about lGeV 2 . To 
achieve this the knowledge of low energy QCD parameters is necessary. 

In order to fix the notations I present here the form of QCD Lagrangian: 

l = i £<(v,7, + im q )r q - ~ g; u g; u , (i) 

q 

where 

V, = d,-ig-Al 

G% = d,Al-d y Al + gf nml A™A l v (2) 

ipq and A™ are quark and gluon fields, a — 1, 2, 3; n, m, I — 1, 2, ...8 are colour indeces, A™ and / are 
Gell-Mann matrices and /-symbols, m q - are bare (current) quark masses, q — u,d, s,c... 

Vacuum condensates are very important in the elucidation of the QCD structure and in description of 
hadron properties at low energies. Condensates, particularly, quark and gluonic ones, were investigated 
starting from the 70-ties. Here, first, it should be noted the QCD sum rule method by Shifman, 
Vainshtein, and Zakharov 1 , which was based on the idea of the leading role of condensates in the 
calculation of masses of the low-lying hadronic states. In the papers of the 70-80-ies it was assumed 
that the perturbative interaction constant is comparatively small (e.g., a s (lGeV) ~ 0.3), so that it 
is enough to restrict oneself by the first-order terms in a s and sometimes even disregard perturbative 
effects in the region of masses larger than 1 GeV. At present it is clear that a s is considerably larger 
(a s (lGeV) ~ 0.5). In a number of cases there appeared the results of perturbative calculations in order 
a 2 and a z s . New, more precise experimental data at low energies had been obtained. 

This review presents the modern status of QCD at low energies. In Chapter 2 the values of light 
quark masses are discussed. Chapter 3 contains the definition of condensates and the description of 
their general properties. In Chapter 4 the QCD coupling constant a s is determined from the data on 
hadronic r-decay and its evolution with Q 2 (in 4-loops approximation) is given. In Chapter 5 quark 
and gluon condensates are found from the r-decay data on V-A, V+A and V correlators. The sum 
rules for nucleon mass with account of a s corrections are analyzed in Chapter 6 and it is shown, that 
they are well satisfyed at the same value of quark condensate, which was found from the r-decay data. 
Various ways of gluon condensate determination: a) from V + A correlators; b) from charmonium sum 
rules are considered in Chapter 7. The QCD sum rules for valence quark distributions in nucleon are 
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presented in Chapter 8, valence u- and c/-quark distributions at low Q 2 were found and the restriction 
on condensates were obtained. Finally, chapter 9 summarizes the state of the art of low energy QCD. 



2 The masses of light quarks 

The u, d, s quark masses had been first estimated by Gasser and Leutwyler about 30 years ago: it was 
demonstrated that m u ,m d ~ 5 MeV and m s ~ lOOMeV PIE]. In 1977 Weinberg 0j, using partial 
conservation of axial current and Dashen theorem j3] to account for electromagnetic selfenergies of 
mesons had proved, that the ratios m u /m d and m s /md may be expressed through K and 7r masses: 

m u _ m 2 K+ - m 2 K0 + 2m^ - 
m d m 2 K0 - m 2 K+ + m\+ 

m s = m 2 K o + m 2 K+ + 
m d m 2 K0 - m 2 K+ + m\+ 

Numerically, (jHJ) and (jlj) are equal 

— !t = 0.56 — = 20.1 5 

m d m d 

Basing on consideration of mass splitting in baryon octet Weinberg assumed, that m s = 150MeV at 
the scale of about 1 GeV. Then 

m u = 4.2MeV, m d = 7.5MeV, m s = lbOMeV (6) 

at 1 GeV. The large m s /m d ratio explains the large mass splitting in pseudoscalar meson octet. For 
m 2 K+ /m 2 r+ we have [fh = (m u + m d )/2] 

rn^ = ms + m = ^ 
m + 2m 

in a perfect agreement with experiment. The ratio m 2 /m 2 expressed in terms of quark mass ratios is 
also in a good agreement with experiment. 

The ratios ©?© were obtained in the first order in quark masses. Therefore, their accuracy is of 
order of accuracy of SU(3) symmetry, i.e. about 20%. 

In jH] it was demonstrated that there is a relation valid in the second order in quark masses 

\m d J Q 2 \m d J 

Using Dashen theorem for electromagnetic selfenergies of it and i^-meson, one may express Q as 

q2 = ( m K« + m K+ ~ m l+ + ™lo)(™K° + m K+ ~ m l+ ~ m lo) / Q N 

D 4ml (m 2 K0 - m 2 K+ + m\+ - m\ Q ) 1 

Numerically, Qd is equal: Qd = 24.2. However, Dashen theorem is valid in the first order in quark 
masses. The electromagnetic mass difference of i^-mesons calculated in [Z| by using Cottingham formula 
and in [S] by large iV c approach increased Am^ = (Mr-+ — Mjfo) e , m , from Dashen value Am^ = 1.27 MeV 
to Arrik = 2.6 MeV and, correspondingly, decreased Q to Q = 22.0 ± 0.6. The other way to find Q is 
from f] — > 7T + 7r~7r° decay, using the chiral effective theory. Unfortunately, the next to leading corrections 
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are large in this approach [§J, what makes uncertain the accuracy of the results. It was found from 
the rj — > 7T + 7r~7r° decay data with the account of interaction in the final state: Q = 22 A ± 0.9 |10j . 
Q = 22.7 ± 0.8 mi and Q = 22.8 ± 0.4 [T2] (the latter from the Dalitz plot). So, the final conclusion is 
that Q is in the interval 21.5 < Q < 23.5. (It must be mentioned that the experiment, where T(i] — > 27) 
and, consequently, T(r] — > 7r + 7r~7r°) were measured by the Primakoff effect, is absent in the last edition 
of the Particle data [T^j, while it persisted in the previous ones. See |T3J for the review.) The ratio 
7 = m u /md can also be found from the ratio of ip' — > (J/ip)rj and ip' — > (J/ip)ii decays [H2 dj . In 
[TK] it was proved that 



F(ip' -»• J/ip + 7T°) „/l-7 



ON /, n 2 / \ 4 / N 3 




(10) 



where and p v are the pion and 77 momenta in ip' rest frame. Eq. (JlO)) is valid in the first order in 
quark mass. The Particle Data Group [T3] gives 

r exp = (3.04 ± 0.71) ■ 10- 2 (11) 

In the recent CLEO Collaboration experiment ^1 ^ was found: r exp = (4.01 ± 0.45) ■ 10~ 2 . Averaging 
these two experimental numbers, assuming the theoretical uncertainty in (jlOj) as 30% and adding in 
quadratures the theoretical and experimental errors, we get from (fTU|) 

Tfl 

7 = — = 0.407 ±0.060 (12) 
m d 

The value close to (|12[) was found recently in I#. The substitution of (fT2"|) into (jHJ) with the account of 
the mentioned above uncertainty of Q, results in 

777- 

— =20.8 ±1.3 (13) 

The value (fT2^1 is slightly lower, then the lowest order result ©, (fTB^l agrees with it. The values (fT ^l .lfTB ^) 
are in agreement with recent lattice calculations [P9*] . 

The calculation of absolute values of quark masses is a more subtle problem. First of all, the masses 
are scale dependent. In perturbation theory their scale dependence is given by the renormalization 
group equation: 

^4 = -7[ a„(A*) f-4 = -f>«V)^ (14) 



In f!14|) a = a s /7T,ji = 1,72 = 91/24, 73 = 10.48 for 3 flavours in MS scheme j20]- I n the first order in 
a s it follows from (pPf]) that: 



m(Q 2 ) 
m(/x 2 ) 



a s (Q 2 



2\ 1 7m 



(15) 



where 7 TO = —4/9 is the quark mass anomalous dimension. There is no good convergence of the series 
(JHJ) below fi = 2 GeV(a s (2 GeV) = 0.31). The recent calculations of m s by QCD sum rules |2*T| . 
from the r-decay data [221 EH] and on lattice (THl 121], are in a not quite good agreement with one 
another. The mean value estimated in is: m s (2 GeV) ~ 105 MeV" with an accuracy of about 
20%. By taking m s (l GeV)/m s (2 GeV) = 1.35 we have then: m s (l GeV) w 142 Me^ and, according 
to m d (l GeV) = 6.8 MeV,m u (l GeF) = 2.8 MeV. The difference m d - m u is equal to: 

n^d — iriu — 4.0 ± 1.0 MeV. This value agrees with one found by QCD sum rules from baryon octet 
mass splitting [25J and D and D* isospin mass differences j^HJj m d ~ m u = 3 ± 1 MeV. For m u ± m d 
we have m u + m d = 9.6 ± 2.5 MeV in comparison with m u = m d = 12.8 ± 2.5 MeV found in [27j. For 
completeness I present here also the value of m c {m c ) (see below, Sec. 7.1) 

m c (m c ) = 1.275 ± 0.015 GeV (16) 
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3 Condensates 



3.1 General properties 

In QCD (or in a more general case, in quantum field theory) by condensates one mean the vacuum mean 
values (0|Oj|0) of the local (i.e. taken at a single point of space-time) of the operators Oj(x), which 
arise due to nonperturbative effects. The latter point is very important and needs clarification. When 
determining vacuum condensates one implies the averaging only over nonperturbative fluctuations. 
If for some operator Oi the non-zero vacuum mean value appears also in the perturbation theory, it 
should not be taken into account in determination of the condensate - in other words, when determining 
condensates the perturbative vacuum mean values should be subtracted in calculation of the vacuum 
averages. One more specification is necessary. The perturbation theory series in QCD are asymptotic 
series. So, vacuum mean operator values may appear due to one or another summing of asymptotic 
series. The vacuum mean values of such kind are commonly to be referred to vacuum condensates. 

In quantum field theory it is assumed, that vacuum correlators Uab(%, y) in coordinate space of any 
two local operators A(x), B(y) 

U AB (x,y) = (0\T{A(x),B(y)}\0) 

at space-like (x — y) 2 < 0, and small x — y (x — y — > 0) may be represented as an operator product 
expansion (OPE) series 

U AB (x-y)=Yl a i( x ~ y)(o\Oi(o)\o), 

i 

where cii(x — y) are called the coefficient functions and are given by perturbation theory. (The strict 
proofs of this statement were obtained only in perturbation theory and for some models). Here, again one 
must take care of separation of perturbative and nonperturbative parts in the definition of condensates. 
The perturbation expansion for aj(x — y) is an asymptotic series and the terms which arise by summing 
of such series may be interpreted as contributions of higher dimension operators. a«(x — y) may be 
infra-red divergent. This is a signal of appearance of an additional condensate in OPE. Also, probably, 
OPE for Hab{x — y) is asymptotic series. In order to avoid all these problems in practical calculations, 
it is necessary to require a good convergence of OPE and perturbation series in the domain of interest. 

Separation of perturbative and nonperturbative contribution into vacuum mean values has some 
arbitrariness. Usually [2E1 EH] , this arbitrariness is avoided by introducing some normalization point 
fi 2 (fi 2 ~ lGeV 2 ). Integration over momenta of virtual quarks and gluons in the region below /i 2 is 
referred to condensates, above fi 2 - to perturbative theory. In such a formulation condensates depend 
on the normalization point fi: (0|Oj|0) = (OlOjlO)^. Other methods for determination of condensates 
are also possible (see below Sec. 5. 2). 

In perturbation theory, there appear corrections to condensates as a series in the coupling constant 

CO 

(0\Oi\0) Q = (OICMO^E C«(Q,/i)<(/x) (17) 

n=0 

The running coupling constant a s at the right-hand part of (fTTj) is normalized at the point /i. The 
left-hand part of (j!7|) represents the value of the condensate normalized at the point Q. Coefficients 
C^(Q,/i) may have logarithms InQ 2 / fi 2 in powers up to n for C®. Summing up of the terms with 
highest powers of logarithms leads to appearance of the so-called anomalous dimension of operators, so 
that in general form it can be written 

<0|O 4 |0) g = E c«(Q, jiKM, (18) 
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where ji - are anomalous dimensions (numbers), and have already no leading logarithms. If there 
exist several operators of the given (canonical) dimension, then their mixing is possible in perturbation 
theory. Then the relations (|17j) . (jl8|) become matrix. 

In their physical properties condensates in QCD have much in common with condensates appearing 
in condensed matter physics: such as superffuid liquid (Bose-condensate) in liquid 4 ife, Cooper pair 
condensate in superconductor, spontaneous magnetization in magnetic etc. That is why, analogously 
to effects in the physics of condensed matter, it can be expected that if one considers QCD at finite 
temperature T, with T increasing at some T = T c there will be phase transition and condensates (or a 
part of them) will be destroyed. Particularly, such a phenomenon must hold for condensates responsible 
for spontaneous symmetry breaking - at T = T c they should vanish and symmetry must be restored. 
(In principle, surely, QCD may have a few phase transitions). 

Condensates in QCD are divided into two types: conserving and violating chirality. As was demon- 
strated in previous Chapter, the masses of light quarks u, d, s in the QCD Lagrangian are small com- 
paring with the characteristic scale of hadronic masses M ~ 1 GeV. In neglecting light quark masses 
the QCD Lagrangian becomes chiral-invariant: left-hand and right-hand (in chirality) light quarks do 
not interact with each other, both vector and axial currents are conserved (except for flavour- singlet 
axial current, non-conservation of which is due to anomaly). The accuracy of light quark masses neglect 
corresponds to the accuracy of isotopical symmetry, i.e. a few per cent in the case of u and d quarks 
and of the accuracy of SU(3) symmetry, i.e. 10-15 % in the case of s-quarks. In the case of condensates 
violating chiral symmetry, perturbative vacuum mean values are proportional to light quark masses and 
are zero within m u = = m s = 0. So, such condensates are determined in the theory much better 
than those conserving chirality and, in principle, may be found experimentally with a higher accuracy. 

Among chiral symmetry violating condensates of the most importance is the quark condensate 
(0|gg|0) (q = u,d are the fields of u and d quarks). (0|gg|0) may be written in the form 

(0\qq\0)^(0\q L q R + q R q L \0) (19) 

where qL,qn are the fields of left-hand and right-hand (in chirality) quarks. As follows from (J19|) . the 
non-zero value of quark condensate means the transition of left-hand quark fields into right-hand ones 
and its not a small value would mean the chiral symmetry violation in QCD. (If chiral symmetry is not 
violated spontaneously, then at small m u ,md (0|gg|0) ~ m U) md). By virtue of isotopical invariance 

(0\uu\0) = (0\dd\0) (20) 

For quark condensate there holds the Gell-Mann-Oakes-Renner relation [30J 

1 m 2 f 2 

(o|ot|o> = -- (21) 

2 m u + m d 

Here m n , f n are the mass and constant of 7r + -meson decay (m n = 140MeV, f n = 131 MeV), m u and 
md are the masses of u and <i-quarks. Relation (}2~T]) is obtained in the first order in m u , m^, m s (for its 
derivation see, e.g. [SI])- To estimate the value of quark condensate one may use the values of quark 
masses m u + = 9.6 MeV, presented in Sec. 2. Substituting these values into (|2Tj) we get 

(0\qq\0) = -(260 MeV) 3 (22) 

The value (JBJ) has characteristic hadronic scale. This shows that chiral symmetry which is fulfilled 
with a good accuracy in the light quark lagrangian (m u ,md/M ~ 0.01), is spontaneously violated on 
hadronic state spectrum. 

An other argument in the favour of spontaneous violation of chiral symmetry in QCD is the existence 
of massive baryons. Indeed, in the chiral-symmetrical theory all fermionic states should be either 
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massless or parity-degenerated. Obviously, baryons, in particular, nucleon do not possess this property. 
It can be shown (321 ETTj . that both these phenomena - the presence of the chiral symmetry violating 
quark condensate and the existence of massive baryons are closely connected with each other. According 
to the Goldstone theorem, the spontaneous symmetry violation leads to appearance of massless particles 
in the physical state spectrum - of Goldstone bosons. In QCD Goldstone bosons can be identified with 
a 7r-meson triplet within m u , m& — > 0, m s ^ (SU(2)-symmetry) or with an octet of pseudoscalar 
mesons (ir, K, rj) within the limit m u , m^, m s —>■ (SU(3)-symmetry). The presence of Goldstone bosons 
in QCD makes it possible to formulate the low-energy chiral effective theory of strong interactions (see 
reviews |SB] , |S1! , [S] ) • 

Quark condensate may be considered as an order parameter in QCD corresponding to spontaneous 
violation of the chiral symmetry. At the temperature of restoration of the chiral symmetry T = T c 
it must vanish. The investigation of the temperature dependence of quark condensate in the chiral 
effective theory shows that (0|gg|0) vanishes at T = T c w 150 — 200MeV. Similar indications were 
obtained also in the lattice calculations [HEj - 

Thus, the quark condensate: 1) has the lowest dimensions (d=3) as compared with other condensates 
in QCD; 2) determines masses of usual (nonstrange) baryons; 3) is the order parameter in the phase 
transition between the phases of violated and restored chiral symmetry. These three facts determine its 
important role in the low-energy hadronic physics. 

Let us estimate the accuracy of numerical value of (|2*2*j) . The quark condensate, as well as quark 
masses depend on the normalization point and have anomalous dimensions equalling to 7^ = — 7 m = |. 
In (|22|) the normalization point fi was taken fi ~ 1 GeV. The Gell-Mann-Oakes-Renner relation is 
derived up to correction terms linear in quark masses. In the chiral effective theory it is possible to 
estimate the correction terms and, thereby, the accuracy of equation (|2*T|) appears to be of order 10%. 
The accuracy of the above taken value m u + = 9.6 MeV which enters (|21j) seems to be of order 
20%. The value of the quark condensate may be also found from the sum rules for proton mass (see 
Chapt.6) as well as from structure functions at r-decay (Chapt.5). The quark condensate of strange 
quarks is somewhat different from (OlwwlO). In |32j it was obtained 

(0|ss|0)/(0|uu|0) = 0.8 ±0.1 (23) 

The next in dimension (d = 5) condensate which violates chiral symmetry is quark gluonic one: 

0<O|^— G£,g|0) =ml(0\qq\0) (24) 

Here - is the gluonic field strength tensor, A n - are the Gell-Mann matrices, a^ v = (i/2)( r y fl r y v — 7^7^). 
The value of the parameter wij was found in [37j from the sum rules for baryonic resonances 

m\ = 0.8 ± 0.2 GeV 2 (25) 

The same value of m^ was found from the analysis of 5-mesons by QCD sum rules [HE], close to (|25j) 
value of TOq = 1.0 GeV 2 was calculated in the model of field correlators [HHI- The anomalous dimension 
of the operator in (|24|) is small [10] ■ Therefore the anomalous dimension of m\ is approximately equal 
to 7 m = -4/9. 

Consider now condensates conserving chirality. Of fundamental role here is the gluonic condensate 
of the lowest dimension: 

(OI^G^JO) (26) 

Since gluonic condensate is proportional to the vacuum mean value of the trace of the energy-momentum 
tensor 9^ u its anomalous dimension is zero. The existence of gluonic condensate had been first indicated 
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by Shifman, Vainshtein, and Zakharov pQ. They had also obtained its numerical value from the sum 
rules for charmonium: 

(0|^G^J0) = 0.012GW 4 (27) 

As was shown by the same authors, the nonzero and positive value of gluonic condensate means, that the 
vacuum energy is negative in QCD: vacuum energy density in QCD is given by e = —(9/32) (0| (a s /7r)G 2 |0) 
The persistence of quark field in vacuum destroys (or suppresses) the condensate. Therefore, if quark 
is embedded into vacuum, this results in its excitation, i.e, in increasing of energy. Thereby, it become 
possible to explain the bag model in QCD: in the domain around quark there appears an excess of 
energy, which is treated as the energy density B in the bag model. (Although, the magnitude of B, 
does not, probably, agree with the value of e which follows from (|27j)). In ref. pQ perturbative effects 
were taken into account only in the order a s , the value for a s being taken about two times smaller 
as the modern one. Later many attempts were made to determine the value of gluonic condensate by 
studying various processes and by applying various methods. But the results of different approaches 
were inconsistent with each other and with ()27j) and sometimes the difference was even very large - the 
values of condensate appeared to be by a few times larger. All of this requires to reanalyse the methods 
of (0|^G 2 |0) determination basing on modern values of a s that will be done in Sections 7,8. 
The d=6 gluonic condensate is of the form 

g 3 f abc (0\G« u G b uX G^\0), (28) 

(f abc - are structure constants of SU(3) group). There are no reliable methods to determine it from 
experimental data. There is only an estimate 0T] which follows from the model of deluted instanton 
gas: 

9 3 f abc (0\G%Gl x G^\0) = ^(12^)1(01^10), (29) 

where p c is the instanton effective radius in the given model (for estimation one may take p c ~ (1/3 — 
l/2)/m). 

The general form of d=6 condensates built from quark fields is: 

a s (0\qiO a qi ■ q k O a q k \0) (30) 

where qi, q k are quark fields of u, d, s quarks, O a - are Dirac and SU (3) matrices. Following [T], Eq.([3*U|) 
is usually factorized: in the sum over intermediate state in all channels (i.e, if necessary, after Fierz- 
transformation) only vacuum state is taken into account. The accuracy of such approximation ~ 1/-^?; 
where N c is the number of colours i.e.~ 10%. After factorization Eq. ()30j) reduces to 

a s (0\m\0) 2 , (31) 

if q = u, d. The anomalous dimension of (|31|) is - 1/9 and it can be approximately put to be zero. And 
finally, d=8 quark condensates assuming factorization reduce to 

as(0|gg|0)-m 2 (0|gg|0) (32) 

(The notation of (}2"^) is used). It should be noted, however, that the factorization procedure in the d=8 
condensate case is not quite certain. For this reason, it is necessary to require their contribution to be 
small. 

There are few gluon and quark-gluon condensates of dimension 8. (The full list of them is given 
in |12]-) As a rule, factorization hypothesis is used for their calculation. The other way to estimate 
the values of these condensate is to use the dilute instanton gas model. However, the latter for some 
condensates gives the results (at accepted values of instanton gas model parameters) by one order of 
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magnitude larger, than the factorization method. The arguments were presented [43J, that instanton 
gas model overestimates the values of d — 8 gluon condensate. Therefore, the estimates based on 
factorization hypothesis are more reliable here. 

The violation of factorization hypothesis is more strong for higher dimension condensates. So, this 
hypothesis may be used only for their estimations by the order of magnitude. 



3.2 Condensates, induced by external fields 

The meaning of such condensates can be easily understood by comparing with analogous phenomena 
in the physics of condensed matter. If the above considered condensates can be compared, for instance 
with ferromagnetics, where magnetization is present even in the absence of external magnetic field, 
condensates induced by external field are similar to dia- or paramagnetics. Consider the case of the 
constant external electromagnetic field F^. In its presence there appears a condensate induced by 
external field (in the linear approximation in F^): 

(0\qa^q\0) F = e gX F^{0\qq\0) (33) 

As was shown in ref.jSj, in a good approximation (Olgcr^glO)^ is proportional to e q - the charge of 
quark q. Induced by the field vacuum expectation value (0\qa^q\0) f violates chiral symmetry. So, it is 
natural to separate (0|gg|0) as a factor in eq. ljHB^l . The universal quark flavour independent quantity % 
is called magnetic susceptibility of quark condensate. Its numerical value had been found in [35] using 
a special sum rule: 

X = -(5.7±0.6)GW 2 (34) 

Another example is external constant axial isovector field A^, the interaction of which with light quarks 
is described by the Lagrangian 

L' = («7a*75^ - dlvlb^A^ (35) 
In the presence of this field there appear induced by it condensates: 

(0\u llllb u\Q) A = -(0\d 7 ^ 5 d\0) A = flA„ (36) 

where f n = 131MeV is the constant of ir — > {iv decay. The right-hand part of eq. ()36|) is obtained 
assuming m u , — > 0, m 2 — > and follows directly from consideration of the polarization operator of 
axial currents Il^(g) in the limit q — > 0, when because of axial current conservation the nonzero con- 
tribution into U^ u (q) q ^o emerges only from one-pion intermediate state. Eq.(|3l)j) was used to calculate 
the axial coupling constant in /3-decay qa [46J. An analogous to ()36|) relation holds in the case of octet 
axial field. Of special interest is the condensate induced by singlet (in flavours) constant axial field 

(0|ji°5 ) |0) = 3/ 2 4°) (37) 

3 { S = «7m75« + 0^75^ + 37^755 (38) 
and the Lagrangian of interaction with external field has the form 

L' = (39) 

Constant fo cannot be calculated by the method used when deriving eq.()36|). since singlet axial current 
is not conserved because of anomaly and the singlet pseudoscalar meson rj' is not Goldstone one. The 
constant /q is proportional to topological susceptibility of vacuum |17j 

fo = ^/V(0), (40) 
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where Nf is the number of light quarks, Nf = 3, and the topological susceptibility of the vacuum x{q 2 ) 
is defined as 

X {q 2 )=il d 4 xe^(0\TQ 5 (x), Q 5 (0)|0) (41) 



QbW = ^Gl u {x)G; u {x\ (42) 

where is dual to G™ : G™^ = (l/2)e jttI/ ,x .G" (T . Using the QCD sum rule, one may relate /g with the 
part of proton spin E, carried by quarks in polarized ep (or ftp) scattering j?7]. The value of /q was 
found from the selfconsistency condition of the obtained sum rule (or from the experimental value of 
S): 

/ 2 = (2.8 ± 0.7) • 10- 2 GeV 2 (43) 

The related to it value of the derivative at q 2 = of vacuum topological susceptibility x'(0), (more 
precisely, its nonperturbative part) is equal to: 

X'(0) = (2.3 ± 0.6) • 10" 3 GeV 2 (44) 

The value is °f essential interest for studying properties of vacuum in QCD. 



4 Test of QCD at low energies on the basis of r-decay data 

4.1 Determination of a s {m 2 r ) 

Collaborations ALEPH |3H] , OPAL and CLEO [HO] had measured with a good accuracy the relative 
probability of hadronic decays of r-lepton R T = B(r — > u T + hadrons) / B{t — > u T eV e ), the vector V and 
axial A spectral functions. Below I present the results of the theoretical analysis of these data basing on 
the operator product expansion (OPE) in QCD ( see a l so jSHlEI!)- In the perturbation theory 

series the terms up to a* will be taken into account, in OPE - the operators up to dimension 8. I 
restrict myself to the case of equal to zero total hadronic strangeness. 
Consider the polarization operator of hadronic currents 

K = * J e^^TJ,(x)J^}dx=(q^-6^q 2 )I^P(q 2 ) + q,q^ ) (q 2 ) } (45) 

where J = V, A; = u^^d, = uj^d. 
The spectral functions measured in r-decay are imaginary parts of Hj(s) and IIj^(s), s = q 2 

vja^s) = 2icImIL§) A (s + iO), a {s) = 2^ImIlf(s + i0) (46) 

Functions Tly\q 2 ) and H A \q 2 ) are analytical functions in the q 2 complex plane with a cut along 
the right-hand semiaxis starting from Am 2 , for Uy\q 2 ) and 9m 2 for IT^^g 2 ). Function U^\q 2 ) has 
kinematical pole at q 2 = 0, since the physical combination, which have no singularities is d^q 2 !^ (q 2 ). 
Because of axial current conservation in the limit of massless quarks this kinematical pole is related to 
one-pion state contribution into LTa^), which has the form [ST] 

f 2 m 2 f 2 

KM* = --^(Wv - <W 2 ) - -fq»q» 2 - 2 (47) 

p qz gz qi _ m z 

The chiral symmetry violation may result in corrections of order f 2 (m 2 /m 2 ) in IT^^g 2 ) (m p is the 
characteristic hadronic mass), i.e. in the theoretical uncertainty in the magnitude of the residue of 
kinematical pole in Il^g 2 ) of order A/ 2 // 2 ~ m 2 /m 2 . 
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Consider first the ratio of the total probability of hadronic decays of r-leptons into states with zero 
strangeness to the probability of r — ► v T ev e . This ratio is given by the equality 

B(t — > v T + hadronss=o) 
Rt,v+a = — 



/ds ( s ^ 
9 1 9 

o T x T 



B(t — > v T ev e 



1 + 2 A, ) (vx + a x + oo)(s) - 2^a (s) 



(4£ 



where \V u d\ = 0.9735±0.0008 is the matrix element of the Kobayashi-Maskawa matrix, Sew — 1.0194± 
0.0040 is the electroweak correction |56j . Only one-pion state is practically contributing to the last term 
in [SB] and it appears to be small: 

AR^ = -247r 2 ^i = -0.008 (49) 

Denote 

u(s) = v 1 + a 1 + a = 27rIm[U { y ] (s) + 1$ (s) + (s)} = 2irImU(s) (50) 

As follows from eq. (|47|) . IL(s) has no kinematical pole, but only right-hand cut. It is convenient to 
transform the integral in eq.([48|) into that over the circle of radius m 2 in the complex s plane [57 -[59]: 

Rr,v+A = G™\V ud \ 2 S EW <f (l + 2^]u(s)+AR^ (51) 

Eq.(|51jl allows one to express R T y + A in terms of n(s) at large | s \— m 2 , where perturbative theory 
and OPE are valid. 

Calculate first the perturbative contribution into eq. (J5T|l . To this end, use the Adler function D{Q 2 ): 
D(Q 2 ) = -2vr 2 ^Q = J2 K na n , « = Q 2 = ~ s , (52) 



the perturbative expansion of which is known up to terms ~ a*. In MS regularization scheme Kq = 
K\ = 1, K'2 — 1.64 [Hill, K% = 6.37 [HI] for 3 flavours and for K± there are the estimates K± = 25 ± 25 
and K4 = 27 ± 16 [63 . The renormgroup equation yields 

= -(3(a) = - £ (3 n a^ 2 (53) 
dlnQ 2 ^ 

1 Q 2 r m da 

/r Ja(p?) f3(a) 



in the MS scheme for three flavours fa = 9/4, (3 X = 4, (3 2 = 10.06, /3 3 = 47.23 [HU E3 EH!- 
Integrating over eq.()52p and using eq.(|33*)l we get 

1 r a (Q 2 ) da 
U(Q 2 ) = — / D(a)^rr (55) 
VV ' 2tt 2 Jay) V ' (3(a) K ' 

Put fi 2 = m 2 and choose some (arbitrary) value a(m 2 ). With the help of eq. (|54j) one may determine 
then a(Q 2 ) for any Q 2 and by analytical continuation for any s in the complex plane. Then, calculating 
(|55|) find n(s) in the whole complex plane. Substitution of n(s) into eq.([51|) determines R T for the 
given a(m 2 ) up to power corrections. Thereby, knowing R T from experiment it is possible to find the 



12 



PT+OPE 


S-plane 

works 






1 ' — 


\ cut 


1 GeV 2 -Q2 


PT+OPE fails 









Figure 1: The applicability region of PT and OPE in the complex plane s. In the dashed region PT + 
OPE does not work. 



corresponding to it a(m 2 ). Note, that with such an approach there is no need to expand the denominator 
in eqs. (J54|) .([55 |) in the inverse powers of lnQ 2 /fi 2 . Advantages of transformation of the integral over 
the real axis (|48|) in the contour integral are the following. It can be expected that the applicability 
region of the theory presented as perturbation theory (PT) + operator product expansion (OPE) in 
the complex s-plane is off the dashed region in Fig.l. It is evident that at positive and comparatively 
small s PT+OPE does not work. 

As is well known, in perturbation theory, in the expansion over the powers of inverse InQ 2 , in the 
first order in 1/lnQ 2 the running coupling constant a s (Q 2 ) has an unphysical pole at some Q 2 = Qq. 
If (3(a) is kept in the denominator in (|54j) . then in ra-loop approximation (n > 1) a branch cut with a 
singularity ~ (1 — Q 2 /Qo)~ 1 ^ n appears instead of pole. The position of the singularity is given by 

<4 = - / #T < 56 > 

/i 2 J (3(a) 

Near the singularity the last term in the expansion of (3(a) dominates and gives the beforementioned 
behavior. Since the singularity became weaker, one may expect a better convergence of series, which 
would allow one to go to lower Q 2 . 

The real and imaginary parts of a s (s)/'K 1 obtained as numerical solutions of eq. (|54jl for various 
numbers of loops are plotted in Fig.2 as functions of s = —Q 2 . The r-lepton mass was chosen as 
normalization point, fi 2 = m 2 and a s (m 2 ) = 0.355 was put in. As is seen from Fig.2, at negative s the 
perturbation theory converges at s < —lGeV 2 and in order to have a good precision of the results 4 
loops calculations are necessary. At positive s, especially for Im(a s /7c), the convergence of the series 
is much better. This comes from the fact, that in the chosen integral form of renormalization group 
equation (|5~4*|) the expansion over n / ln(Q 2 / A 2 ) is avoided, this expansion being not a small parameter 
at intermediate Q 2 . (The systematical method of analytical continuation from the spacelike to timelike 
region with summation of tt 2 terms was suggested in \§T\ and developed in jHEl)- F° r instance, in the 
next to leading order 



2nImU(s + iO) = 1 ' 



ir(3 
instead of 



7T 

— — arctg 



(57) 



2nlmll(s + i0) = l+ (58) 

(3 ln(s/A 2 ) 

which would follow in the case of small n/ln(s/A 2 ). 

The a s (Q 2 ) at Q 2 > in low Q 2 region (0.8 < Q 2 < 5GeV 2 ) is plotted in Fig.3. (Four loops 
are accounted, a s (m 2 ) is put to be equal to a s (m 2 ) = 0.33. As follows from r-decay rate a s (m 2 ) = 
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Figure 2: Real and imaginary parts of cx^(s)/7r as an exact numerical solution of RG equation 
on real axes for different number of loops. The initial condition is chosen a s = 0.355 at s = — m 2 ., 
Nf = 3. Vertical dotted lines display the position of the unphysical singularity at s = —Qq for each 
approximation (4 — > 1 from left to right). 



0.352 ± 0.020 and the value of one standard deviation below the mean one is favoured by low energy 
sum rules). 

Integration over the contour allows one to obviate the dashed region in Fig.l (except for the vicinity 
of the positive semiaxis, the contribution of which is suppressed by the factor (1 — ) 2 in eq.([5T|)). i.e. 
to work in the applicability region of PT+OPE. 

The OPE terms, i.e., power corrections to polarization operator, are given by the formula PQ: 



nonpert 



(0 2n ) 



71 



+ 



M • - 

n>2 

128 , . _ , , 2 
ixa s {0\qq\ 0) 



81Q 



7_as 

6 7T 

'29 17, Q 



Q 4 



(m u + m d ) (0 | qq | 0) 



a. 



1 + — + rr: hi — 



24 18 yi 



+ 



Q 8 



(59) 



(a s -corrections to the 1-st and 3-d terms in eq. (f5T)|) were calculated in (HHl and [10], respectively). Con- 
tributions of terms proportional to m 2 , m 2 d are neglected. When calculating the d=6 term, factorization 
hypothesis was used. Gluon condensate of dimension d = 6 g 3 (0 \ G 3 \ 0) ^Hj) does not contribute 
to polarization operator ()59|) . This is a consequence of the general theorem, proved by Dubovikov and 
Smilga [7T], that in case of self-dual gluonic fields there are no contributions of gluon condensates of 
dimensions higher than d = 4 to vector and axial currents polarization operators. Since the vacuum 
expectation value of G 3 operator does not vanish for self-dual gluonic fields, this means the vanishing 
of the coefficient in front of g 3 (0 \ G 3 \ 0) condensate in (|5T?|) . The same argument refers to dimension 8 
gluon operators g A G A with the exception of some of them, like g A [ G^Gfy — (l/4)5 a/ 5G^G™ 1/ ] 2 , which 
have zero expectation values in any self-dual field. But the latter are suppressed by a small factor 
1/47T 2 arising from loop integration in comparison with tree diagram, corresponding to d — 8 four quark 
condensate (0$) ~ (qGq ■ qq) contribution. The contribution from this condensate may be estimated 
as | (0$) \< 10~ 3 GeV |S2] (see below, Sec. 5.1) and appears to be negligibly small. The d = 8 two 
quarks - two gluons operator 0' 8 ~ g 2 DqGGq is nonfactorizable, its vacuum mean value is suppressed 
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Figure 3: a s (Q 2 ) normalized at m^, a s {m 2 ) = 0.33. 

by 1/N C and one may believe, that its contribution to (|59|) is also small. It can be readily seen that d=4 
condensates (up to small a s corrections) give no contribution into the integral over contour eq. (|51j) . 
Rt,v+a may be represented as 

Rr.y+A = S\V ud \ 2 S EW ( 1 + 5' em + 5 (0) + 5${ A ) + Ai? {0) = 3.486 ± 0.016 (60) 

where 5' em = (5/l27r)a em (m 2 ) = 0.001 is the electromagnetic correction [72] . S^]_ v = — (3.3 ± 1.1) ■ 10~ 3 
is the contribution of d=6 condensate (see below) and 6^ is the PT correction. The right-hand part 
presents the experimental value obtained as a difference between the total probability of hadronic decays 
R T = 3.647 ± 0.014 [73] and the probability of decays in states with the strangeness S = — 1 R TjS = 
0.161 ± 0.007 [HI EE]- For perturbative correction it follows from eq. (l60j) . that 

5 {0) = 0.208 ±0.006 (61) 

From (joTj) employing the above described method, the constant a s {m 2 T ) was found [S2] 

a s (m 2 T ) = 0.352 ±0.020 (62) 

The calculation was made with the account of terms ~ a^, the theoretical error was assumed to be equal 
to the last term contribution. May be, the error is underestimated (by ~ 0.010), since the theoretical 
and experimental errors were added in quadratures. The value a s {m 2 ) (jo^j) corresponds to: 

a s (m 2 z ) = 0.121 ±0.002 (63) 

This value is in agreement with recent determination [7^] of a s {m 2 ) from the whole set data 

a s (m 2 z ) = 0.1182 ±0.0027 (64) 

4.2 Instanton corrections 

Some nonperturbative features of QCD may be described in the so called instanton gas model (see [77] 
for extensive review and the collection of related papers in |7J|). Namely, one computes the correlators 
in the S'?7(2)-instanton field embedded in the SU (3) color group. In particular, the 2-point correlator 
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of the vector currents had been computed long ago jHDj- Apart from the usual tree- level correlator 
~ In Q 2 it has a correction which depends on the instanton position and radius p. In the instanton gas 
model these parameters are integrated out. The radius is averaged over some concentration n(p), for 
which one or another model is used. Concerning the 2-point correlator of charged axial currents, the 
only difference from the vector case is that the term with 0-modes must be taken with opposite sign. In 
coordinate representation the answer can be expressed in terms of elementary functions, see jHU]. An 
attempt to compare the instanton correlators with the ALEPH data in the coordinate space, was made 
inRef.[EU. 

We shall work in momentum space. Here the instanton correction to the spin- J parts Iv> J ' of the 
correlator (}4l)j) can be written in the following form: 



1. illSt(? 



r(0) 

1 A inst 



dpn(p) 
dpn(p) 



4 

Sq 1 
4 

7 ' 



I 2 2 


1/2 y 


[-p q 


0,0, -2 | 



n d) 



n (!) (J\ rr(0) 



V, inst 







(65) 



Here K\ is modified Bessel function, G^(z\ . . .) is Meijer function. Definitions, properties and approx- 
imations of Meijer functions can be found, for instance, in [H2]- In particular the function in ()65j) can 
be written as the following series: 



30 
13 



1/2 ' 
0,0,-2 



4 

3? 



2 i « k r(fc + i/2) 
z 20F ^ r 2 (fc + i)r(fc + 3) 



x 



{ [lnz + i[)(k + 1/2) - 2i/)(k + l) - ij(k + 3)] 2 
+ ip'(k + 1/2) - 2 ip'(k + 1) - i)'(k + 3) } (66) 

where ip{z) = T'(z)/T(z). For large \z\ one can obtain its approximation by the saddle-point method: 



^13 



1/2 
0,0,-2 



^-3/2 e -2Vi 



\z\ > 1 



(67) 



The formulae should be treated in the following way. One adds n inst to usual polarization 
operator with perturbative and OPE terms. But the terms ~ 1/g 4 must be absorbed by the operator 
4 in Eq.(fE5|). since the gluonic condensate (G 2 ) is averaged over all field configurations, including 
the instanton one. Notice negative sign before 1/g 4 in Eq. (J64|) . This happens because the negative 
contribution of the quark condensate (mqq) in the instanton field exceeds positive contribution of the 
gluonic condensate (G 2 ). In real world (mqq) is negligible. 

The correlators ([53} possess appropriate analytical properties, they have a cut along positive real 
axes: 



Imn« inst (g 2 + *0) 
(0) (g 2 + *0) = 



1/2 
0,0,-2 



imn 



A, inst \ 



dpn(p) 



2tt 2 P 2 



We shall consider below the instanton gas model. It is a model with fixed instanton radius 

n(p) = n 5(p - p ) 

In it was estimated: 



(68) 
(69) 

(70) 
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1/3 fm w 1.5 - 2.0 GeV" 



no 



1 fm 



;i.O - 1.5) x 10- J GeV 



(71) 
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Figure 4: The instanton correction to the r decay ratio versus po (a) and "versus r mass" (b) for no = 
1.5 x 10~ 3 GeV 4 . The thin solid lines in FigEjD are the values of 1+6^ (sq), where 6^(sq) are perturbative 
corrections, calculated as described in Sec. 4.1. The upper curve corresponds to a s {m 2 ) = 0.355, the 
lower one - to a s (m^) = 0.330. The shadowed regions represent the uncertainties in perturbative 
calculations, the dark shadowed band is their overlop. The dashed lines are 1 + <r°)(s) + $inst,8 (s) 
corresponds to a s (m 2 ) = 0.330. 



In fact, the instanton liquid model, with the account of instanton self- interaction was mainly considered 
in [77j, but the arguments, from which the estimations (J71|) follow, refer also to the instanton gas model. 
In this case, the value of n (f7T|) should be considered as an upper limit (see also f7Sj). 

Now we consider the instanton contribution to the r-decay branching ratio. Since the instanton 
correlator ()65|) has 1/q 2 singular term in the expansion near (see Eq. (|66|) ) . the integrals must be 
taken over the circle, as in (joTj) . In the instanton model the function ao(s) differs from experimental 5- 
function, which gives small correction. So we shall ignore the last term in (J48)) and consider the integral 
with IIy^ A + in (Jo"Tj) . The instanton correction to the r-decay branching ratio can be brought to 
the following form: 

5inst = -48^ r d pn(p)p i GH Uml 1/2 ) « sin(2p m T ) (72) 

jo y u, — i, — 4 j Po m T 

Since the parameters (J7TJ) are determined quite approximately, we may explore the dependence of <5i nst 
on them. The <5i nst versus po for fixed no = 1.5 ■ !0~ 3 GeV 4 is shown in Fig0] 

As seen from FigEK the instanton correction to hadronic r-decay is extremely small except for 
unreliably low value of the instanton radius po < 1.5 GeV -1 . At the favorable value [ZZj po = 1.7 GeV -1 
the instanton correction to R T is almost exactly zero. This fact confirms the calculations of a s (m 2 ) 
(Sec. 4.1), where the instanton corrections were not taken into account. 

Eq. ()72|) can be used in another way. Namely, the r mass can be considered as free parameter so- The 
dependence of the fractional corrections 5^ and ^olLo + ^mst on s is shown in FigE^. The result strongly 
depends on the instanton radius and rather essentially on the density n . For p = 1.7 GeV -1 and 
no = lfm -4 , the instanton curve is outside the errors already at So ~ 2 GeV 2 , where the perturbation 
theory is expected to work. 

We came to the conclusion, that in case of variable r mass the instanton contribution becomes large 
at so < 2GeV 2 . That means, that -R Ti v + a(so) given by (foTj) cannot be represented by PT+OPE at 
so < 2GeV 2 and the results, obtained in this way are not reliable. 
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4.3 Comparison with other approaches 

There are many calculations of a s (m 2 ) from the total r-decay rate, using the same idea, which was 
used above - the countour improved fixed order perturbation theory [531 EHl EZJ-jSHl EH]- (For more 
recent ones, see |23[ I54|.) The results of these calculations coincide with presented above in the limit of 
errors and give a s (m 2 ) = 0.33 — 0.35. From these values by using renormalization group one can find 
a s {m 2 ) = 0.118 — 0.121 in agreement with a s {m 2 z ) determinations from other processes (see [IB], [IB])- 

Till now only one renormalization scheme was considered - the MS scheme. In BLM renormalization 
scheme [83J, which have some advantages from the point of view of perturbative pomeron theory [84J, 
the result is a s (m 2 ) = 0.621 ± 0.008 [EH], corresponding in the framework of BLM scheme to the same 
value of a s {m 2 ) = 0.117— 0.122. At low scales, however, the a s (Q 2 ) behavior is essentially different 
from that, presented in Fig. 3. 

Few words about a s calculations in analytical QCD (see [HE] and references herein). According to 
this theory the coupling constant a s (Q 2 ) is calculated by renormalization group in the spacelike region 
Q 2 > 0. Then, by analytical continuation to s = —Q 2 > Ima s (s) was found on the right semiaxes. 
It was assumed, that a s (s) is an analytical function in the complex s-plane with a cut along the right 
semiaxes < s < oo. The analytical a s (s) an is then defined in the whole s-plane by dispersion relation. 
Such a s (s) an has no unphysical singularities. Let us calculate a s {m 2 ) an using the same experimental 
data as before, i.e. 5^ given by Eq.(|61|). In the analytical QCD the countour integral (JoTj) is equal to 
the integral of ImU(s) over real positive axes. (In the previous calculation the integral was running 
from s = —Qq to m 2 .) Qualitatively, it leads to much smaller R T in the analytical QCD than in the 
conventional approach with the same a s (m 2 ), or vice versa, it is necessary to have much larger a s (m 2 ) an 
in order to get experimental R T . The calculation of integral (joTj) with H(s) expressed through a s (s) ana [, 
shows that experimental R T results in a s (m 2 ) = 0.141 ±0.004 in contradiction with other data. (In [S7] 
an attempt was made to get an agreement of analytical QCD with common value of a s (m 2 ). For this 
goal the constituent quark model with specific quark- ant iquark potential was used in the domain of low 
and intermediate s. Evidently, such approach cannot be considered as a s determination in QCD: in 
this approach QCD is modified on large circle in complex plane of the radius | s \— m 2 in contradiction 
with the basic assumption of a s calculation from hadronic r-decay rate.) 



5 Determination of condensates from spectral functions of r- 
decay 

5.1 Determination of quark condensate from V — A spectral function 

In order to determine the quark condensate from r-decay data it is convenient to consider the difference 
V — A of polarization operators Ily^ — , where the contribution of perturbative terms is absent. 
]Iy^(s) — n^(s) is represented by OPE: 

n<!» W -n« W = E^i(i + CD ^)) (73) 

The gluonic condensates contribution drops out in the V — A difference and only the following conden- 
sates up to D=10 remain 

0\~ A = 2 K + m d ) (0 | qq \ 0) = - f 2 m 2 n [1] (74) 
0\- A = 2ira s (0 | (u lfl X a d)(d lfl X a u) - (u l5lfl X a d)(d l5lfl X a u) | o) = 

18 



(75) 



0\~ A = -8Tca s m 2 (0 | qq | 0) 2 , jHHl EH 1 (76) 

(77) 



OX, A = -^a s (0 | qq | 0) 2 



^ + 32tt 2 (0|^G 2 |0) 

9 71 



where m 2 , is determined in eq.(J2U). In the right-hand of (|75|) . (|75j) . (|77jl the factorization hypothesis 
was used. For 6 operator it is expected pQ, that the accuracy of factorization hypothesis is of order 
1/iV 2 ~ 10%, where N c = 3 is the number of colours. For operators of dimensions d > 8 the factor- 
ization procedure is not unique. (But, as a rule, the arising differences are not very large - for d — 8 
operator entering eq.fJZBJ) it is about 20%). The accuracy of factorization hypothesis becomes worse 
with increasing of operator dimensions: for Og~ A , it is worse, than for 0\~ A and for OY Q ~ A it is worse 
than for OX~ A . 

Operators O4 and 0§ have approximately zero anomalous dimensions, the Og anomalous dimension 
is equal to - 11/27. Calculations of the coefficients in front of a s in eq. (J7H|) gave C4 = 4/3 (HOI and 
c 6 = 89/48 [91]. (For A the a 2 correction is known [20]: (59/6)(a s /vr) 2 ).) The a s corrections to 0%~ A 
are unknown - they are included into the not certainly known value of m 2 ,, a s corrections to Oiq are 
unknown also. (In this Section indeces V — A will be omitted and Or, will mean condensates with a s 
corrections included.) 

Our aim is to compare the OPE theoretical predictions with the experimental data on V —A structure 
functions measured in r-decay and with the help of such comparison to determine the magnitude of the 
most important condensate Oq. The condensate O4 is small and is known with a good accuracy: 

4 = -0.5 ■ 10" 3 GeV 4 (78) 

We put uIq = 0.8Ge^ 2 and in the analysis of the data the values of the condensates 0$ and Oiq are 
taken to be equal to 

8 = -2.8 • 10~ 3 GeV 8 (79) 
Oio = -2.6 • 10" 3 GeV 10 (80) 

and their Q 2 -dependence, arising from anomalous dimensions is neglected. 

In the calculation of numerical values (|75 j) . (|79*|) it was assumed, that aq q (l GeV 2 ) = 
= -(2tt) 2 (0 I qq | 0) 1GeV = 0.65 GeV 3 , (0 | (a s /n)G 2 | 0) = 0.005 GeV 4 - see below, eq.'s flSTjl . (fTTTj) . 

As was shown in jHIj the dimension d = 8 four-quark operators for vector and axial currents are 
of opposite sign and equal in absolute values up to terms of order 1/N 2 : Og = —0 A (1 + 0(N~ 2 )). 
(The exact value of N~ 2 correction is uncertain - it depends on factorization procedure.) So, for 0^ +A 
we have from (|79*j) the estimation: | 0% +A \< 10~ 3 GeV 8 , which was used in calculation U(s) nonpert , 
Eq.((HSD- 

For riy^(s) — n^(s) substractionless dispersion relation is valid: 

U^(s) - ri?(s) = JL r Vl{t) - ai{t) dt + & m) 



Z7T Z JO t 



(The last term in the right-hand part is the kinematic pole contribution). The experimental data for 
Vi(s) — ai(s) are presented in Fig. 5 

In order to improve the convergence of OPE series as well as to suppress the contribution of large s 
domain in dispersion integral use the Borel transformation. Put s = s^e 1 ^ (<ft = on the upper edge of 



1 There was a sign error in the contribution of 0$ in |5T] , 
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Figure 5: The measured difference v\{s) — a\{s). Figures from [IB] and |%9] . reproduced in [5T] . 

the cut) and make the Borel transformation in sq. As a result, we get the following sum rules for the 



real and imaginary parts of (J81|) : 
s 



s \ ds 
( OS| sin0J — ai)(s) — 



k=l 



k cos (k(p)0 



2k+2 



k\M 2k 



exp 



M 2 



COS( 



s \ ds 



E(-i) 

k=l 



fc sin(/c0) 2k +2 



k\ M 2AH 



(82) 



(83) 



The use of the Borel transformation along the rays in the complex plane has a number of advantages. 
The exponent index is negative at ir/2 < <p < 3tt/2. Choose in the region tt/2 < <p < tt. In this region, 
on one hand, the shadowed area in Fig. 1 in the integrals is touched to a less degree, and on 

the other hand, the contribution of large s, particularly, s > m 2 , where experimental data are absent, 
is exponentially suppressed. At definite values of <fi the contribution of some condensates vanishes, what 
may be also used. In particular, the condensate 0$ does not contribute to ([H2~)) at (f> = 5rr/6 and to 
(I83|) at <p = 27r/3, while the contribution of Oq to (J82j) vanishes at cf) = 3tt/A. Finally, a well known 
advantage of the Borel sum rules is factorial suppression of higher dimension terms of OPE. Figs. 6, 7 
present the results of the calculations of the left-hand parts of eqs.(JB2])?(JEIIl) on the basis of the ALEPH 
[15] experimental data comparing with OPE predictions - the right-hand part of these equations. 

When comparing the theoretical curves with experimental data it must be taken in mind, that the 
value of which in the figures was taken to be equal to experimental one f v = 130.7MeV , in fact has 
a theoretical uncertainty of the order (Af 2 / f 2 ) theor ~ m 2 /m 2 , where m p is characteristic hadronic scale 
(say, p-meson mass). This uncertainty is caused by chiral symmetry violation in QCD. Particularly, the 
account of this uncertainty may lead to a better agreement of theoretical curve with the data in Fig.6b. 
The calculation of instanton contributions fEq. ()65[) L shows, that in all considered above cases they are 
less than 0.5 ■ 10~ 3 at M 2 > 0.8 GeV 2 , i.e. are well below the errors. (In some cases they improve the 
agreement with the data.) The best fit of the data (the dashed curves at Fig.'s 6,7) was achieved at 
the value 

6 = -4.4 • 10~ 3 GeV 6 (84) 
It follows from ([54^1 after separating a s correction [1 + (89/48)a s /7r) = 1.33]: 



a s (Q | qq | 0) 2 = 1.5 • 10~ 4 GeV t 



J5) 
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Figure 6: Eq. (|82j) : the left-hand part is obtained basing on the experimental data, the shaded region 
corresponds to experimental errors; the right-hand part - the theoretical one - is represented by the 
dotted curve, numerical values of condensates 4 , 8 , O w , 6 are taken according to (|75|) . (|75jl . (|5nj) . (|5Hl : 

a) = 3vr/4, b) = 5tt/6. 



The error may be estimated as 30%. The value (|K5|) in the limit of errors agrees with previous estimation 
[51j . The contribution of dimension 10 is negligible in all cases at M 2 > 1 GeV 2 . It is worth mentioning 
that the theory, i.e. the OPE agrees with the data at M 2 > 0.8 GeV 2 . The good agreement of 
the theoretical curves with the data confirms the chosen value of Os ()78)1 and, therefore, the use of 
factorization hypothesis. From (84), with the use of a s (lGeV 2 ) = 0.55 (see Fig.3) the value of quark 
condensate at 1 GeV can be found 

(0|gg|0)i GeV = -1-65 ■ 10~ 2 GeV 3 = -(254 MeVf (86) 

and the convenient parameter is 

a 9q (l GeV 2 ) = -(27r) 2 (0|gg|0) 1 GeV = 0.65 GeV 3 (87) 

The magnitude of quark condensate (86) is close to that which follows from the Gell-Mann-Oakes-Renner 
relation (eq.flHJ)). 

In the last years there were many attempts [2I|>[EI2]-|IIE| to determine quark condensates using V-A 
spectral functions measured in r-decay. Unlike the approach presented above, where the polarization 
operator analytical properties were exploited in the whole complex g 2 -plane, what allowed one to sep- 
arate the contribution of operators of different dimensions, the authors of jEB],[H2]-|nE] considered the 
finite energy sum rules - FESR (or integrals over contours) with chosen weight functions. In [021 HH] 
the N c — > oo limit was used. In |9*3*1 |03 103 EH] an attempt was made to find higher dimension conden- 
sates (up to 18 in [87], up to 16 in [013 El and up to 12 in [95]). Determination of higher dimension 
condensates requires fine tunning of the upper limit of integration in FESR. If the upper limit of inte- 
gration so in FESR is below 2 GeV 2 (e.g., such an upper limit, so = 1.47 GeV was chosen in j0B|), then 
instanton-like corrections, not given by OPE are of importance. (See Sec. 4. 2). The same remark refers 
to the case of weight factors singular at s = 0, like s~ l , I > [53:, when there is an enhancement of the 
contribution of low s, where OPE breaks down. Taking in mind these remarks, we have a satisfactory 
agreement of the values of condensate (jMJ), presented above, with those found in [SHI EH EE] • 
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Figure 7: The same for eq. (J83|) : a) (p = 2tt/3, b) (f> = 3n / 4. 



5.2 Determination of condensates from V + A and V structure functions of t -decay 

Let us turn now to study the V + A correlator in the domain of low Q 2 , where the OPE terms play a 
much more essential role, than in the determination of R T . A general remark is in order here. As was 
discussed in Ref.|2*H] and stressed recently by Shifman the condensates cannot be defined in rigorous 
way, because there is some arbitrariness in the separation of their contributions from perturbative part. 
Usually [23 EH] they are defined by introduction of some normalization point /j, 2 with the magnitude of 
few Aq Cd . The integration over momenta in the domain below /i 2 is addressed to condensates, above /i 2 - 
to perturbation theory. In such formulation the condensates are /i-dependent (On) = (Od)^ and, strictly 
speaking, they also depend on the way how the infrared cut-off /i 2 is introduced. The problem becomes 
more severe when the perturbative expansion is performed up to higher order terms and the calculation 
pretends on high precision. Mention, that this remark does not refer to chirality violating condensates, 
because perturbative terms do not contribute to chirality violating structures in the limit of massless 
quarks. For this reason, in principle, chirality violating condensates, e.g. (0|gg|0), can be determined 
with higher precision, than chirality conserving ones. Here I use the definition of condensates, which 
can be called n-loop condensates. As was formulated in Chapt.4, we treat the renormalization group 
equation (J54j) and the equation for polarization operator (J55j) in n-loop approximation as exact ones; 
the expansion in inverse logarithms is not performed. Specific values of condensates are referred to such 
procedure. Of course, their numerical values depend on the accounted number of loops; that is why the 
condensates, defined in this way, are called n-loop condensates. 

Consider the polarization operator II = riy+ A + IT^, defined in (I5U|) and its imaginary part 

uj(s) = Ui(s) + ai(s) + a (s) = 27rImII(s + zO) (88) 

In parton model uj(s) — * 1 at s —>■ oo. Any sum rule can be written in the following form: 

J f(s) uj cxp (s) ds = Z7r j f(s) n thcor (s) ds (89) 

where f(s) is some analytical in the integration region function. In what follows we use o; exp (s), obtained 
from r-decay invariant mass spectra published in j^H] for < s < m 2 . The experimental error of the 
integral is computed as the double integral with the covariance matrix uj(s)uj(s') —ZJ(s)uJ(s'), which 
also can be obtained from the data available in Ref . j3H] • in the theoretical integral in (|H9*j) the contour 
goes from so + iO to So — iO counterclockwise around all poles and cuts of theoretical correlator U(s). 
Because of Cauchy theorem the unphysical cut must be inside the integration contour. 
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Figure 8: The results of the Borel transformation of V + A correlator for two values a s (m 2 ) = 0.355 
and a s (m 2 ) = 0.330. The widths of the bands correspond to PT errors, dots with dashed errors - 
experimental data. The dashed curve is the sum of the perturbative contribution at a s (m 2 ) = 0.330 
and 4 , (Eq.'s (59),fJ21)) and 6 (Eq.'s ©,©) condensate contributions. 



The choice of the function f(s) in Eq.()89|) is actually a matter of taste. At first let us consider usual 
Borel transformation: 

B„ P (M 2 ) = [' e-/"^ p M |L = B pt (M>) + ^ Y. ( J";;^„ PO) 

We separated out the purely perturbative contribution B pt , which is computed numerically according 
to (J89j) and Eqs. (j52|) - (J55|) . Remind that Borel transformation improves the convergence of OPE series 
because of the factors l/(n — 1)! in front of operators and suppresses the contribution of high-energy 
tail, where the experimental error is large. But it does not suppress the unphysical perturbative cut, 
the main source of the error in this approach, even increase it since e~ s l M > 1 for s < 0. So the 
perturbative part B pt (M 2 ) can be reliably calculated only for M 2 > 0.8 — 1 GeV 2 and higher; below 
this value the influence of the unphysical cut is out of control. 

Both B exp and B pt in 4-loop approximation for a s (m 2 ) = 0.355 and 0.330 are shown in Fig. 8. The 
shaded areas display the theoretical error. They are taken equal to the contribution of the last term in 
the perturbative Adler function expansion K^a 4 (p)2*|) . 

The contribution of the 0$ operator is of order Og~ A /N 2 and negligible [5T]- (In fact, it depends 
on the factorization procedure and uncertain for this reason). The contributions of D = 4 and D = 6 
operators are positive [see So, the theoretical perturbative curve must go below the experimental 

points. The result shown in Fig. 8 is in the favour of the lower value of the QCD coupling constant 
a s (m 2 ) = 0.330 (or, may be, a s {m 2 ) = 0.340). As is seen from Fig.8, the theoretical curve (perturbative 
at a s (m 2 ) = 0.330 plus OPE terms) is in agreement with experiment at M 2 > 0.9 GeV 2 . 

In order to separate the contribution of gluon condensate let us perform the Borel transformation 
along the rays in the complex s-plane in the same way, as it was done in Sec. 5.1. The real part of the 
Borel transform at <fi = 5ir/6 does not contain d = 6 operator. 

ReB exp (M 2 e 1 ^ 6 ) = ReB pt (M 2 e^ 6 ) + n 2 ^ (91) 

The results are shown in Fig.9. If we accept the lower value of a s (m 2 ), we get the following restriction 
on the value of gluon condensate: 

G a ^ u G a ^ = 0.006 ± 0.012 GeV 4 , a s (m 2 T ) = 0.330 and M 2 > 0.8 GeV 2 (92) 
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Figure 9: Real part of the Borel transform (91) along the ray at the angle (f) = 5n/6 to the real axes. 
The dashed line corresponds to the gluonic condensate given by the central value of (|92J1 . 

The theoretical and experimental errors are added in quadratures in Eq. (j92)l . 

Turn now to analysis of the vector correlator (the vector spectral function was published by ALEPH 
in jin]). In principle this cannot give any new information in comparison with V — A and V + A 
cases. However the analysis of the vector current correlator is important since it can also be performed 
with the experimental data on e + e~ annihilation. The imaginary part of the electromagnetic current 
correlator, measured there, is related to the charged current correlator (|45jl by the isotopic symmetry. 
The statistical error in e + e~ experiments is less than in r-decays because of significantly larger number 
of events. So it would be interesting to perform similar analysis with e + e~ data, which is a matter for 
separate research. 

At first we consider usual Borel transformation for vector current correlator, since it was originally 
applied in |l()()j for the sum rule analysis. It is defined as (89) with the experimental spectral function 
^ex P = 2di instead of v\ + a% + ao (the normalization is v\(s) — > 1/2 at s — ► oo in parton model). 
Respectively, in the r.h.s. one should take the vector operators 20 v = O v+A + O v ~ A . The numerical 
results are shown in Fig. 10. The perturbative theoretical curves are the same as in Fig.8 with V + A 
correlator. The dashed lines display the contributions of the gluonic condensate given by Eq. (j92J) . 
201 = - 3 - 5 x 10" 3 GeV 6 and 20% = 0%~ A = -2.8 x 10~ 3 GeV 8 added to the a s (m 2 T ) = 0.330- 
perturbative curve. The contribution of each condensate is shown in the box below. Notice, that for 
such condensate values the total OPE contribution is small, since positive O4 compensate negative Oq 
and Og. The agreement is observed for M 2 > 0.8 GeV 2 . 

The Borel transformations along the rays in the complex plane results in the same conclusion; at 
M 2 > 0.8 — 0.9 GeV 2 the agreement with experiment at 2% level is achieved at a s (m 2 ) = 0.33 — 0.34 
and at the values of quark and gluon condensates given by (84) and There is some discrepancy in 
the vector spectral function obtained in r-decay and in e + e~ annihilation (see [51] . |101j and references 
herein): the e + e~ data are below r-decay ones by 5-10% in the interval s = 0.6 — 0.8 GeV 2 . The 
substitution of e + e~ data instead of r-decay data in the sum rule presented in Fig. 10 does not spoil the 
agreement of the theory with experiment in the limit of errors. 

A few words about instanton contributions. They can be calculated in the same way, as in the 
case of V — A correlators. At the chosen values of instanton gas parameters instanton contributions 
are small, less than 0.5 • 10~ 3 at M 2 > 0.8 GeV 2 and do not spoil the agreement of the theory with 
experiment. 
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Figure 10: Borel transformation for vector currents. 



6 Determination of quark condensate from QCD sum rules 
for baryon masses 

Since in QCD with massless quarks the baryon masses arise due to spontaneous violation of chiral 
symmetry and in a good approximation, the proton mass (as well as A-isobar) can be expressed through 
quark condensate [32], the QCD sum rules for baryon masses are a suitable tool for determination of 
quark condensate assuming that baryon masses are known. The sum rules can be derived by considering 
the polarization operator 

n(p)=i / d 4 xe ipx (0\T{r]{x),fj{0)}\0} (93) 



where r](x) is the quark current with baryon quantum numbers. In case of proton the most suitable 
current is |3~2] 11772]. 

■nix) = e abc (u a C^u b ) l5 ^d c (94) 

where u a , d c - are u and d quark fields, a, b, c are colour indeces, C is the charge conjugation matrix. 
After Borel transformation the sum rules for proton mass have the form [321 EH IS] 



M 6 E 2 (s /M 2 )L- 4/9 
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Here M is the Borel parameter, m is the nucleon mass, aq q is given by (j87j). 7^ = 0.577. 
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Figure 11: The sum rules for proton mass Eq.s' (fHKJ) . (jnflj) . The dashed and dash-dotted curves give 
A 2 , determined correspondingly from (J9"5j) and ()96|) . the experimental value of m was substituted (left 
scale). The solid line gives m as the ratio of (}9l)j) to (jHSJ). 



(99) 



L - corresponds to anomalous dimensions, Sq is the continuum threshold and fi 2 is the normalization 
point, chosen as fi 2 = 1 GeV 2 . The constant X p is defined as A 2 = 2(2-7r) 4 A 2 



(0\r)\p) 



Apfp, 



(100) 



where v p is the proton spinor. The a s corrections to proton sum rules were found in |103j . The function 
f(s) is small, |/| < 0.2 at 0.9 < M 2 < 1.5 GeV 2 and a s correction to the term proportional to a? g 
can be neglected. The sum rules (JHHJ), (J9l)|) were calculated at the following values of parameters: 
(0\(a s /n)G 2 \0) = 0.005 GeV 4 , {b = 0.20 GeV 4 ), s = 2.5 GeV 2 . The numerical value of quark 
condensate was not fixed by the value given in (|H7|) . but considered as a free parameter. For the best 
fit of the sum rules it was chosen to be aq g = 0.60 (cf.fjHZJ)). First, the values of X 2 was found from 
(l93 j) . (pm j) . where the experimental value of proton mass was substituted, Fig. 11, left scale. Then Eq. (j9l)j) 
was divided by (|93j) and the theoretical value of the proton mass was found, Fig. 11, right scale. 

As is seen from Fig.ll, A 2 , determined from ()95|) . (|9l)|) are almost independent on M 2 and coincide 
with one another, as it should be. The proton mass value coincide with the experimental one with a 
precision better than 3%. The conclusion is, that the value 



"qq 



0.60 GeV 



(101) 



describes well the proton mass sum rule. The main source of the error is the large a s correction (about 
0.8) to the first term in (j9*3j) . If we suppose that its uncertainty is 20%, then the corresponding error 

Therefore, we get from proton mass sum rules 



in aq q is ±0.1 GeV 



-qq 



(0.60 ±0.10) GeV 3 



(102) 



A remark about a possible role of instantons in the sum rules for proton mass. As was found in 
|. |105j if the quark current with proton quantum numbers is given by (J9"5j) . then instantons do not 
change the sum rule Their contribution to (96) is moderate in instanton gas model, if the model 
parameters are chosen as in (fTTj) [10 ^1105) and may shift the value of quark condensate (102) by 10-20%, 
i.e. in the limit of quoted error. 
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7 Gluon condensate and determination of charmed quark mass 
from charmonium spectrum 



7.1 The method of moments. The results 

The existence of gluon condensate had been first demonstrated by Shifman, Vainshtein and Zakharov 
PP. They considered the polarization operator IT c (g 2 ) of the vector charmed current 

n c (9 2 )(^-W) = * / d 4 xe iqx (0\TJ^x), J v (0)\0) (103) 

J M (x) = C7^c (104) 

and calculated the moments of II c (g 2 ) 

An 2 d 

M n (Q 2 ) = ^r(-^) n n c (Q 2 ), (105) 

(Q 2 = —q 2 ) at Q 2 = 0. The OPE for H(Q 2 ) was used and only one term in OPE series was accounted - 
the gluonic condensate. In perturbative part of H(Q 2 ) only the first order term in a s was accounted and a 
small value of a s was chosen, a s (m c ) ~ 0.2. The moments were saturated by contribution of charmonium 
states and in this way the value of gluon condensate (J27)) was found. The SVZ approach pQ was criticized 
in |1061| . where it was shown that the higher order terms of OPE, namely, the contributions of G 3 and G A 
operators are of importance at Q 2 = 0. Reinders, Rubinstein and Yazaki |107j demonstrated, however, 
that SVZ results may be restored, if one considers not small values Q 2 > instead of Q 2 = 0. Later 
there were many attempts to determine the gluon condensate by considering various processes within 
various approaches. In some of them the value (J27j) (or ones, by a factor of 1.5 higher) was confirmed 
|1001 11081 HO^j , in others it was claimed, that the actual value of the gluon condensate is by a factor 
2-5 higher than flUD [HU] . 

From today's point of view the calculations performed in pQ have a serious drawback. Only the first 
order (NLO) perturbative correction was accounted in PQ and it was taken rather low value of a s , later 
not confirmed by the experimental data. The contribution of the next, dimension 6, operator G 3 was 
neglected, so the convergence of the operator product expansion was not tested. 

There are recent publications |lllj where the charmonium as well as bottomonium sum rules were 
analyzed at Q 2 = with the account of a 2 perturbative corrections in order to extract the charm and 
bottom quark masses in various schemes. The condensate is usually taken to be or some another fixed 
value. However, the charm mass and the condensate values are entangled in the sum rules. This can be 
easily understood for large Q 2 , where the mass and condensate corrections to the polarization operator 
behave as some series in negative powers of Q 2 , and one may eliminate the condensate contribution to 
a great extent by slightly changing the quark mass. Vice versa, different condensate values may vary 
the charm quark mass within few per cents. (See Fig. 12 below.) 

Therefore, in order to perform reliable calculation of gluon condensate by studying the moments 
of charmed current polarization operator it is necessary to account a 2 perturbative corrections to the 
moments, a s corrections to gluon condensate contribution, (G 3 ) term in OPE and to find the region in 
(n,Q 2 ) space, where all these corrections are small. This program was realized in Ref. |112| . The basic 
points of this consideration are presented below. 

The dispersion representation for II(g 2 ) has the form 

R(s) = 47rImII c ( S + *0) , U c (q 2 ) = t- f° R (g) ds , (106) 

Air 2 J Ami s \ s — Q ) 
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where R(oo) = 1 in partonic model. In approximation of infinitely narrow widths of resonances R(s) 
can be written as a sum of contributions from resonances and continuum 

3 7T 

R ^ = r>2 2 { \ J2 m ^ T ^^ S ( s - m %) + 0(s-s o ) (107) 

where Q c = 2/3 is the charge of charmed quarks, So - is the continuum threshold (in what follows 
■s/sq = 4.6 GeV), a(s) - is the running electromagnetic constant, a(m^) = 1/133.6. The polarization 
operator moments are expressed through R as: 

According to (jl(J8|) the experimental values of moments are determined by the equality 

bW j ~ 4«l K + Q 2 ) n+1 "(*o + Q 2 ) n { U9j 

In the sum in (jl(J9|) the following resonances were accounted: J/ip(lS), tp(2S), -0(3770), ?/>(4040), 
•0(4160), ■0(4415), their T^ ee widths were taken from PDG data [13 . It is reasonable to consider 
the ratios of moments M n i(Q 2 ) / M n 2(Q 2 ) from which the uncertainty due to error in Tj/^^ ee markedly 
falls out. Theoretical value for II(g 2 ) is represented as a sum of perturbative and nonperturbative 
contributions. It is convenient to express the perturbative contribution through R(s), making use of 

dinnD , dinHi) : 

R( s ) = J2R (n) {s,fi 2 )a n {^) (110) 

n>0 

where a(/i 2 ) = a s (/i 2 )/7r. Nowadays, three terms of expansion in (jllO}) are known: R<® [EE! r{1) dH, 



R^ [115j . They are represented as functions of quark velocity v — yl — 4m 2 /s, where m c - is the 
pole mass of quark. Since they are cumbersome, I will not present them here (see |112j for details). 

Nonperturbative contributions into polarization operator have the form (restricted by d=6 opera- 
tors): 

IW rt (Q 2 ) = j^y 2 (0 I ~G 2 I 0)[ fi%) + afW(z) ]+ 

+ j^9 3 f abc (0\G; u G b uX Gl II \0)F(z), z = --^ 2 (in) 

Functions / (0) (z), f {1) (z) and F(z) were calculated in pQ, |116j . |117j . respectively. The use of 
the quark pole mass is, however, inacceptable. The matter is that in this case the PT corrections to 
moments are very large in the region of interest and perturbative series seems to diverge. 

So, it is reasonable to use MS mass m(/i 2 ), taken at the point fi 2 = m 2 . The calculations, performed 
in ref.100 show, that in the region near the diagonal in (Q 2 ,n) plane, Q 2 /4m 2 = n/5 — 1 all mentioned 
above corrections are small. For example, 



n = 10,Q 2 = 4m 2 : m = 0.045 , ^ = 1-136 , = -1.673 (112) 

(here mean the coefficients at the contributions of terms ~ a k to the moments, M^ G ' k ' - are the 
similar coefficients for gluonic condensate contribution). 

At a ~ 0.1 and at the ratios of moments given by ()112|) there is a good reason to believe that the 
PT series well converges. Such a good convergence holds (at n > 5) only in the case of large enough 
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Figure 12: The dependence of m(m) on (0|a s /7r)G 2 |0) obtained at n = 10, Q 2 = 0.98 • 4m 2 and 
a s (Q 2 + m 2 ). 

Q 2 , at Q 2 = one does not succeed in finding such n, that perturbative corrections to the moments, a s 
corrections to gluonic condensates and the term ~ (G 3 ) contribution would be simultaneously small. 

It is also necessary to choose the scale - normalization point fi 2 where a s (fi 2 ) is taken. In ()110|) R(s) 
is a physical value and cannot depend on fi 2 . Since, however, we take into account in ()110|) only three 
terms, at unsuitable choice of /i 2 such fj, 2 dependence may arise due to neglected terms. At large Q 2 the 
natural choice is fi 2 = Q 2 . It can be thought that at Q 2 = the reasonable scale is ii 2 = m 2 , though 
some numerical factor is not excluded in this equality. That is why it is reasonable to take interpolation 
form 

fi 2 = Q 2 + m 2 , (113) 

but to check the dependence of final results on a possible factor at m 2 . Equalling theoretical value of 
some moment at fixed Q 2 (in the region where M^> and are small) to its experimental value one 
can find the dependence of m on ((a s /n)G 2 ) (neglecting the terms ~ (G 3 )). Such a dependence for 
n = 10 and Q 2 /Am 2 = 0.98 is presented in Fig. 12. 

To fix both m and ((a s /ir)G 2 ) one should, except for moments, take their ratios. Fig. 13 shows the 
value of m obtained from the moment M 10 and the ratio M 10 /M 12 at Q 2 = 4m 2 and from the moment 
M 15 and the ratio M 15 /M 17 at Q 2 = 8m 2 . The best values of masses of charmed quark and gluonic 
condensate are obtained from fig. 13: 

m(m 2 ) = 1.275 ±0.015 GeV , {~ G2 ) = 0.009 ± 0.007 GeV 4 (114) 

The calculation shows, that the influence of continuum - the last term in eq. ()107p is completely neg- 
ligible. Up to now the corrections ~ (G 3 ) were not taken into account. It appears that in the region 
of n and Q 2 used to find m and gluonic condensate they are comparatively small and, practically, not 
changing m, increase ((a s /ii)G 2 ) by 10 — 20% if the term ~ (G 3 ) is estimated according to (f2~H|) at 
p c = 0.5/m. 

It should be noted that improvement of the accuracy of r j/^-, ee would make it possible to precise the 
value of gluonic condensate: the widths of horizontal bands in fig. 13 are determined mainly just by this 
error. In particular, this, perhaps, would allow one to exclude the zero value of gluonic condensate, that 
would be extremely important. Unfortunately, eq.(114) does not allow one to do it for sure. Diminution 
of theoretical errors which determine the width of vertical bands seems to be less real. 



1 1 1 i 1 1 1 i 1 1 1 1 1 1 i 1 1 1 i 1 1 1 
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Figure 13: The dependence of m(m) on (0|(a s /7r)G 2 |0) obtained from the moments (horizontal bands) 
and their ratios (vertical bands) at different a s . The left-hand figure: Q 2 = 4m 2 , n = 10, M10/M12; 
the right-hand figure - Q 2 = 8m 2 , n = 15, M 15 /M 17 . 

In order to check the results (114) for gluon condensate the pseudoscalar and axial- vector channels 
in charmonia were considered. The same method of moments was used and the regions in the space 
(n, Q 2 ) were found, where higher order perturbative and OPE terms are small. In the pseudoscalar 
case it was obtained |118j that, if for m the value (114) is accepted and the contribution of (0|G 4 |0) 
condensate may be neglected, then there follows the upper limit for gluon condensate 



<0|— G 2 |0) < 0.008 GeV 4 (115) 

71 

The contribution of D = 6 condensate (0|G 3 |0) is shown to be small. If (G 4 ) condensate is accounted and 
its value is estimated by factorization hypothesis, then the upper limit for gluon condensate increases 
to 

<0|— G 2 |0) < 0.015 GeV 4 (116) 

71 

In |119j the case of the axial- vector channel in charmonia was investigated and very strong limitations 
on gluon condensate were found: 

(0|^ G 2 |0) = 0.005±g-g8i GeV 4 (117) 

71 

Unfortunately, ()117|) does not allow one to exclude the zero value for gluon condensate. It should be 
mentioned, that the allowed region in (n, Q 2 ) space, where all corrections are small, is very narrow in 
this case, what does allow us in |119j to check the result ()117j) by studying some other regions in (n, Q 2 ), 
as it was done in the two previous cases - vector and pseudoscalar. 

Let us now turn the problem around and try to predict the width Tj/^ ee theoretically. In order to 
avoid the wrong circle argumentation we do not use the condensate value just obtained, but take the 
limitation (^G 2 ) = 0.006 ± 0.012 GeV 4 found from r-decay data. Then, the mass limits m = 1.28 — 
1.33 GeV can be found from the moment ratios exhibited above, which do not depend on V j/^^ ee if the 
contributions of higher resonances is approximated by continuum (the accuracy of such approximation 
is about 3%). The substitution of these values of m into the moments gives 

rSSUe = 4.9±0.8keV (118) 
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in comparison with experimental value Tj/^ ee = 5.26 ± 0.37 keV. Such a good coincidence of the 
theoretical prediction with experimental data is a very impressive demonstration of the QCD sum rules 
effectiveness. It must be stressed, that while obtaining (|118|) no additional input were used besides the 
condensate restriction taken from Eq.fjHU) and the value of a s {m 2 ). 



7.2 The attempts to sum up the Coulomb-like corrections. Recent publications 

Sometimes when considering the heavy quarkonia sum rules the Coulomb-like corrections are summed 
up [120 -[124]. The basic argumentation for such summation is that at Q 2 = and high n only small 
quark velocities v < 1/ y/n are essential and the problem becomes nonrelativistic. So it is possible to 
perform the summation with the help of well known formulae of nonrelativistic quantum mechanics for 
|-?/>(0)| 2 in case of Coulomb interaction (see pL25j). 

This method was not used here for the following reasons: 

1. The basic idea of our approach is to calculate the moments of the polarization operator in QCD 
by applying the perturbation theory and OPE (l.h.s. of the sum rules) and to compare it with the 
r.h.s. of the sum rules, represented by the contribution of charmonium states (mainly by J/ip in vector 
channel). Therefore it is assumed, that the theoretical side of the sum rule is dual to experimental one, 
i.e. the same domains of coordinate and momentum spaces are of importance at both sides. But the 
charmonium states (particularly, J/ip) are by no means the Coulomb systems. A particular argument 
in favor of this statement is the ratio Y j/^ ee /Y^i^ ee = 2.4. If charmonia were nonrelativistic Coulomb 
system, T^ ee would be proportional to (^(O)) 2 ~ l/(n r + l) 3 , and since ip' is the first radial excitation 
with n r = 1, this ratio would be equal to 8 (see also |125j ). 

2. The heavy quark- ant iquark Coulomb interaction at large distances r > r con f ~ 1 GeV -1 is screened 
by gluon and light quark- ant iquark clouds, resulting in string formation. Therefore the summation of 
Coulombic series makes sense only when the Coulomb radius rc ou i is below r con f. (It must be taken 
in mind, that higher order terms in Coulombic series represent the contributions of large distances, 
r ^ ^Coui-) For charmonia we have 

rcoui « — I « 4 GeV" 1 (119) 

m c C F a s 

It is clear, that the necessary condition -Rcoui < -Rconf is badly violated for charmonia. This means that 
the summation of the Coulomb series in case of charmonium would be a wrong step. 

3. The analysis is performed at Q 2 /4rh 2 > 1. At large Q 2 the Coulomb corrections are suppressed 
in comparison with Q 2 = 0. It is easy to estimate the characteristic values of the quark velocities. 



At large n they are v ~ y(l + Q 2 /4m 2 ) /n. In the region (n,Q 2 ) the exploited above quark velocity 

v ~ 1/ \fb ~ 0.45 is not small and not in the nonrelativistic domain, where the Coulomb corrections are 
large and legitimate. 

Nevertheless let us look on the expression of R c , obtained after summation of the Coulomb corrections 
in the nonrelativistic theory [126J. It reads (to go from QED to QCD one has to replace a — > Cf®. Si 
C 1/3): 

„ 3 7rCVa s 3 / x x 2 x 4 \ , 

Rc Coul = — = -v 1 + - + + ... 120 

ccoui 2 1-e- x 2 V 2 12 720 / K J 

where x = ~nCpa s jv. At v = 0.45 and a s ~ 0.26 the first 3 terms in the expansion ()120j) . accounted in 
our calculations, reproduce the exact value of -R Cj coui with the accuracy 1.6%. Such deviation leads to 
the error of the mass fa of order (1 — 2) x 10 -3 GeV, which is completely negligible. In order to avoid 
misunderstanding, it must be mentioned, that the value of -R c coui, computed by summing the Coulomb 
corrections in nonrelativistic theory has not too much in common with the real physical situation. 
Numerically, at chosen values of the parameters, -R c Coui ~ 1.8, while the real value (both experimental 
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and in the perturbative QCD) is about 1.1. The goal of the arguments, presented above, was to 
demonstrate, that even in the case of Coulombic system our approach would have a good accuracy of 
calculation. 

At v — 0.45 the momentum transfer from quark to antiquark is Ap ~ 1 GeV. (This is a typical 
domain for QCD sum rule validity.) In coordinate space it corresponds to Ar 9 g ~ 1 GeV -1 . Comparison 
with potential models [126J demonstrates, that in this region the effective potential strongly differs from 
Coulombic one. 

4. Large compensation of various terms in the expression for the moments in MS scheme is not 
achieved, if only the Coulomb terms are taken into account. This means, that the terms of non- 
Coulombic origin are more important here, than Coulombic ones. 

For all these reasons the summation of nonrelativistic Coulomb corrections is inadequate in the 
problem in view: it will not improve the accuracy of calculations, but would be misleading. 

In the recent publication |127j it is claimed, that gluon condensate is much larger than the presented 
above values, it was found (0|(a s /7r) G 2 |0) = 0.062 ± 0.019 GeV 4 . The author of |127j considered the 
model, where hadronic spectrum is represented by infinite number of vector mesons. The polarization 
operator, calculated in this model was equalled to those in QCD, given by perturbative and OPE terms. 
The value of gluon condensate was found from this equality. The zero width approximation was used 
for vector mesons. It is clear, however, that the account of non-zero widths results in the terms of the 
same type, proportional to 1/Q 4 , as the contribution of gluon condensate. The sign of these terms is 
such, that they lead to diminishing of gluon condensate. Namely, after accounting for p-meson width, 
the value of gluon condensate decreases by a factor of 2. For this reason the results of [127J are not 
reliable. 



8 Valence quark distributions in nucleon at low Q 2 and the 
condensates 

Quark and gluon distributions in hadrons are not fully understood in QCD. QCD predicts the 
evolution of these distributions with Q 2 in accord with the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) 128J- 130 equations, but not the initial values from which this evolution starts. The stan- 
dard way of determination of quark and gluon distributions in nucleon is the following 131 -[135J (for 
the recent review see |136p . At some Q 2 = Qq (usually, at low or intermediate Q 2 ~ 2 — 5GeV 2 ) the 
form of quark (valence and sea) and gluon distributions is assumed and characterized by the number 
of free parameters. Then, by using DGLAP equations, quark and gluon distributions are calculated at 
all Q 2 and x and compared with the whole set of the data on deep inelastic lepton-nucleon scattering 
(sometimes also with prompt photon production, jets at high p_|_ etc). The best fit for the parameters 
is found and, therefore, quark and gluon distributions are determined at all Q 2 , including their initial 
values q(Qo,x), g(Ql,x). Evidently, such an approach is not completely satisfactory from theoretical 
point of view - it would be desirable to determine the initial distribution directly from QCD. In QCD 
calculation valence quark distributions in nucleon essentially depend on vacuum condensate, particu- 
larly, on gluon condensate. Therefore, the comparison of valence quark distributions calculated in QCD 
with those , found by the fit to the data, allows one to check the values of condensates obtained by 
consideration of quite different physical phenomena. For all these reasons it is desirable to find quark 
and gluon distribution in nucleon at low Q 2 ~ 2 — 5GeV 2 basing directly on QCD. 

The method of calculation of valence quark distributions at low Q 2 (Q 2 = 2 — 5 GeV 2 ) was sug- 
gested in |137j and developed in [138 -[140J. Recently, the method had been improved and valence 
quark distributions in pion |141j and transversally and longitudinally polarized p-meson |142j had been 
calculated, what was impossible in the initial version of the method. The idea of the approach (in the im- 
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proved version) is to consider the imaginary part (in s-channel) of a four-point correlator U(pi,p2, q, q') 
corresponding to the non-forward scattering of two quark currents, one of which has the quantum 
numbers of hadron of interest (in our case - of proton) and the other is electromagnetic (or weak). 
It is supposed that virtualities of the photon q 2 ,q' 2 and hadron currents p\,p\ are large and negative 
I? 2 1 = \l' 2 \ ^ bill \p\\ ^ -R^ 2 ? where R c is the confinement radius. It was shown in (137j that in this 
case the imaginary part in s-channel [s = (pi + q) 2 } of U(pi,p2', q, q 1 ) is dominated by a small distance 
contribution at intermediate x. (The standard notation is used: x is the Bjorken scaling variable, 
x = —q 2 j1v , v = piq). The proof of this statement is given in ref. |138| . So, in the mentioned above 
domain of q 2 , q 12 , pi, pi and intermediate x Imli{p, p 2 ; q, q') can be calculated using the perturbation 
theory and the operator product expansion in both sets of variables q 2 = q' 2 and pl,p\- Only the lowest 
twist terms, corresponding to condition | p\jq 2 |<C 1, | p\jq 2 |<^C 1 are considered. 

The approach is inapplicable at small x and x close to 1. This can be easily understood for physical 
reasons. In deep inelastic scattering at large \q 2 \ the main interaction region in space-time is the light- 
cone domain and longitudinal distances along the light-cone are proportional to 1/x and become large 
at small x |143t 1144] . For OPE validity it is necessary for these longitudinal distances along light-cone 
to be also small, that is not the case at small x. At 1 — x <C 1 another condition of applicability of the 
method is violated. The total energy square s = Q 2 {l/x — 1) +p\ Q 2 = —q 2 is not large at 1 — x -C 1- 
Numerically, the typical values to be used below are Q 2 ~ 5 GeV 2 , pi ~ — 1 GeV 2 . Then, even at 
x ~ 0.7, s ~ 1 GeV 2 , i.e., at such x we are in the resonance, but not in the scaling region. So, one 
may expect beforehand, that our method could work only up to x ~ 0.7. The inapplicability of the 
method at small and large x manifests itself in the blow-up of higher order terms of OPE. More precise 
limits on the applicability domain in x will be found from the magnitude of these terms. 

The further procedure is common for QCD sum rules. On one hand the four-point correlator 
Xl{pi,p 2 ; q, q') is calculated by perturbation theory and OPE. On the other hand, the double dispersion 
representation in pi, pi in terms of physical states contributions is written for the same correlator and 
the contribution of the lowest state is extracted using the Borel transformation. By equalling these two 
expression the desired quark distribution is found. Valence quark distributions in proton according to 
this method were calculated in |145j . The basic results of |145j are presented below. 

Consider the 4-current correlator which corresponds to the virtual photon scattering on the quark 
current with quantum number of proton: 



{p\,p2,q,q) = — i J d A xd 4 yd A z ■ e 



■(0\T{ V (x), j u /(y),j u A0), fj(z)}\0), (121) 

where 7](x) is the three-quark current ()94|) . Choose the currents in the form ft = u^^u, = d'y^d, 
i.e. as an electromagnetic current which interacts only with u(d) quark (with unit charges). Such a 
choice allows us to get sum rules separately for distribution functions of u and d quarks. Let us take the 
hadronic currents momenta to be nonequal, perform the independent Borel transformation over pi and 
p\ and only at very end put the Borel parameters M\ and Mf to be equal. The described procedure 
allows one to kill nondiagonal transitions matrix elements of the type 

(0\j h \h*}(h*\j; l (y)fj(0)\h}(h\j h \0} (122) 
and thus makes it possible to separate the diagonal transition of interest 

(0\j h \h}(h\j; l (y)fj(0)\h}(h\j h \0). (123) 
As was shown in Ref. 11451 the sum rules for nucleon have the form 



2 \2 

JL_ xq u > d (x)e- m2/M2 = ImT° d + Power corrections (124) 



4M 4 32tt 4 
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Figure 14: Bare loop diagrams, corresponding to unit operator contribution for u- and <i-quarks (re- 
spectively, a) and b)). 

Here the l.h.s is the phenomenological side of the sum rule - the proton state contribution, X p is defined 
in (100). In numerical calculations it will be put equal X 2 = 3.0 GeV 6 . The right hand side is calculated 
in QCD. The excited states contribution - the continuum is identifyed with the contribution of bare 
loop diagram, starting from continuum threshold value sq and is transferred to the l.h.s. of the sum 
rule. The bare loop contribution to the sum rules is represented in Fig. 14. 
The results after the double Borel transformation are 



ImTl (d) 



u(d) t s M 2 
^° V ; 32vr 



E 2 (s /M 2 



where 



<p%(x) = x{l - x) 2 {l + 8x), ip d Q {x) = x(l - xf{l + 2x), 
E 2 (z) is given by (98) The substitution of eq.(125) into the sum rules (124) results in 

(d) 2M 6 e m2 / M2 u(d) ( s \ 

X( l( x V = J2 \x) ■ E 2 \j^ J 



(125) 
(126) 

(127) 



Making use of relation X 2 e m2 / M2 = M 6 E 2 which follows from the sum rule for the nucleon mass (see 
(95)) in the same approximation), we get 



i i 
do(x)dx =1, J uo(x)dx = 2 

ii o 



128) 



In the bare loop approximation there also appears the sum rule for the second moment: 



x(q%(x) + qo(x))dx = 1 



(129) 



Analogously to |138j one can show that relations (|128|) . (jl29|) hold also when taking into account 
power corrections proportional to the quark condensate square in the sum rules for the 4-point correlator 
and in the sum rules for the nucleon mass. Relations ()128|) reflect the fact that proton has two u- 
quarks and one <i-quark. Relation ()129|) expresses the momentum conservation law - in the bare 
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loop approximation all momentum is carried by valence quarks. Therefore, the sum rules 
demonstrate that the zero order approximation is reasonable. 

Let us calculate the perturbative corrections to bare loop and restrict ourselves by the leading order 
(LO) corrections proportional to lnQ\j hi 2 , where Ql is the point, where the quark distributions q(x, Qq) 
is calculated and /z 2 is the normalization point. In our case it is reasonable to choose /i 2 to be equal to 
the Borel parameter /z 2 = M 2 . The results take the form: 

d LO 



IX) 



d {x)\l + \ln{Ql/M 2 ).^r 
3 Zn 



l/2 + x + ln((l-x) 2 /x) + 



-5 - I7x + 16x 2 + 12x 3 (3 - 2x)x 2 ln(l /x) 



u LO (x) 



Uq{X 



1 + 



6(1 -x)(l + 2x) 
4 a s (Q 2 



3 2tt 
7 - 59x + 46x 2 + 48x 3 



HQl/M 2 ) 
(15-- 



(1 -x) 2 (l + 2x) 
1/2 + x + ln(l -x) 2 /x+ 
x)x 2 ln(l/x) 



(130) 



6(1 -x)(l + 8 



x 



X 



) 2 (1 + 8x) 



(131) 



where u (x) and d (x) are bare loop contributions, given by (|127|) . 

The contributions of gluon condensate to u and <i-quarks distribution were found to be (the ratios 
to bare loop contributions are presented): 



u(x) {G 2) _ ((a s /n)G 2 ) ii 2 (ll + 4x-31x 2 ) 



u {x) 
d (x) 



M 4 12 x(l -x) 2 {l + 8x) 
((a s /7t)G 2 )7t 2 (l-2x 2 ) 



M 4 



6 x 2 {l- x) 2 (l + 2x) 



E (s /M 2 )/E 2 (^- 2 ] 



E (s /M 2 )/E 2 (^- 2 ) 



(132) 



(133) 



The contributions of quark condensate - the terms, proportional to a s (M 2 )(0 | qq \ 0) 2 are few times 
smaller, than the contributions of gluon condensate and are not presented here. (They can be found in 
Ref.[HS]). 

The final result for valence quark distribution in proton are of the form 



XU[X) 



M 6 e m 2 /M 2 
A 7V 



2x{l 



x) 2 (l 



8x)Eo 



So . 

M 2 ' 



1 + 



u 



LO, 



x, Q 



u [x) 



+ 



xd(x) 



M 6 e m 2 /M 2 



A 2 



N 



u {x) 
2x(l - 
1 



[u(x) {G 2 } +u(x) aa{qg) 2 



X 



s 



[l + 2x)E 2 (^ 



d (x) 



M 2 ' I 

[d(x) (G 2) + d{x) as( g q) 2] 



1 + 



d LO {x,Ql) 



dn(x) 



+ 



(134) 



(135) 

()135|) for various ((a s /n)G 2 ) = 



The valence u and d quark distribution calculated according to 
0.00,0.006,0.012 GeV A are shown in Fig.15. 

The following values of parameters were used: a s (0 \ qq \ 0) 2 given by (jHB]) . A 2 = 3.0 GeV 6 , 
Aqcd = 250 GeV, so = 2.5 GeV 2 , Ql = 5 GeV 2 . The contribution of (G 3 ) condensate was estimated 
on the basis of instanton model - Eq. (}2*9~j) . This contribution may influence u and o?-quark distributions 
at x < 0.2 and increase both of them by 10-20%. The limits of applicability of QCD calculations 
are: for u- quark - 0.2 < x < 0.65, for ci-quark - 0.3 < x < 0.65. The lower limit arises from gluon 
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Figure 15: d- and w-quark distribution at various values of gluon condensate (((a s /n)G 2 ) = 0.012,0.06 
and GeV 4 , respectively dotted, solid and dashed lines). Thick solid line corresponds to the results of 

usu. 

condensate contribution - it was required, that this contribution does not exceed 30%, the upper limit 
was determined by increasing of perturbative corrections and ct s (qq) 2 terms. For comparison, Fig. 15 
presents the results of the fit to data basing on the solution of DGLAP equation. The LO order fit j!31j 
is chosen, but not the more precise NLO fits 132 -[136J, because QCD calculations of quark distributions 
were performed in LO. As is seen from Fig. 15 the calculated in QCD initial (at Qq = 5 GeV 2 ) valence 
quark distributions at 0.3 < x < 0.7 are in a satisfactory agreement with those found from the data. 
The preferred value of gluonic condensate is (0 | (a s /n)G 2 | 0) = 0.006 GeV 4 , the values higher than 
0.012 GeV 2 and lower than 0, probably, may be excluded. These statements are in full accord with 
those obtained in the previous Chapters. 



9 Conclusion 

The basic parameters of QCD - a s (Q 2 ) and the values of vacuum condensates, determining hadron 
physics at low momentum transfers (Q 2 ~ 1 — 5 GeV 2 ), are reliably determined by the theory. It was 
demonstrated, that the values of these parameters, found by consideration of various processes are in 
a good agreement with one another. The values of u, d, s quark masses and their ratios are known 
now with a good precision - about 10-15% in the ratios and about 20% in the mass absolute values. 
The precision in charmed quark mass value m c (m c ) in MS renormalization scheme is extremingly high 
- about 1%. The knowledge of a s (Q 2 ) and condensates makes it possible to find the polarization 
operators of vector and axial currents at Q 2 > 1 GeV 2 with high precision. In such calculation high 
order perturbative terms - ~ a 2 and, in some cases, ~ oP s must be accounted. Therefore, we have now a 
good basis for theoretical description of many physical phenomena in low energy QCD - hadron masses, 
their static properties, quark distributions in hadrons etc. Of course, at even lower momentum transfer, 
Q 2 < 1 GeV 2 , the approach, exploited in this review and based on perturbation theory and OPE, 
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does not work: the confinement mechanism and the mechanism of spontaneous symmetry breaking are 
acting in full strength. The construction of various models is inavoidable here. But for such models the 
knowledge of basic QCD parameters is also quite important - they may play the role of cornerstones 
for the models. 

I summarize here the final values of a, and condensates: 



a s (m 2 T ) = 0.340 ±0.015 MS - scheme (136) 



a 



(0 I —G* | 0) = 0.005 ± 0.004 GeV 1 (137) 



7T 



(0 I qq | 0)i GeV = -(1.65 ±0.15) ■ 10~ 2 GeV 3 , q = u,d (138) 

a s (0 | qq | 0) 2 = (1.5 ± 0.2) • 10" 4 GeV 6 (139) 

(In determination of ()139j) the factorization hypothesis is assumed.) The values of errors, given in 
f ()137|) - (jl39j) are a bit uncertain, since the procedure of averaging of errors in different processes is 
subjective in essential way. 
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